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ABSTRACT 

Direct imaging of extrasolar planets with future space-based coronagraphic telescopes may provide 
a means of detecting companion moons at wavelengths where the moon outshines the planet. We 
propose a detection strategy based on the positional variation of the center of light with wavelength, 
“spectroastrometry.” This new application of this technique could be used to detect an exomoon, to 
determine the exomoon’s orbit and the mass of the host exoplanet, and to disentangle of the spectra 
of the planet and moon. We consider two model systems, for which we discuss the requirements 
for detection of exomoons around nearby stars. We simulate the characterization of an Earth-Moon 
analog system with spectroastrometry, showing that the orbit, the planet mass, and the spectra of both 
bodies can be recovered. To enable the detection and characterization of exomoons we recommend that 
coronagraphic telescopes should extend in wavelength coverage to 3 micron, and should be designed 
with spectroastrometric requirements in mind. 

Subject headings: planets and satellites: detection — astrometry — techniques: imaging spectroscopy 


1. INTRODUCTION 

Techniques and ideas for detecting moons orbiting exo¬ 
planets have progressed rapidly in the last decade. These 
advances are driven, at least in part, by a desire to ex¬ 
tend our understanding of satellite formation and origins 
to cases beyond the Solar System. For example, while 
the formation of moons and/or moon systems may be 
common around giant planets, there exist many physical 
models for this process jLunine fc Stevenson|1982| Canup 

fc Ward|2QQ2||Mc)squeira ^ Estrada|2UU3[|Ogihara &: ida 

2Q12[). Additiona lly, the likelihoods of satell ite formation 

by giant im pact dHartm ann fc Davis] 1975 
Ward 1976) and captureJ ARCor^iybol 


m 


_ Cameron fc 

j^innon|l984|) 


remain uncerta in (e.g.,|Agnor Asphaug|2UU4||Agnor 

Hamilton||2QQ6 ) 


Exomoons can be potentially habitable targets in their 
own right, and may also affect the habitability an d char¬ 
acterization of the host planet (Heller et al. 12014). Hab- 

)Ia 


itable. Earth-like moons orbiting gas giant plan ets in cir 
cumstellar habitable zones have been proposed ( Williams 
|1997l |Kaltenegger| \2000\j 

lUeller'fc Bai 

tF - 


et al. 


an d discussed dScharf 


Barnes||2Q13| [Heller ^ Fu- 


dritz 2Q15a| |Heller fc Barnes|2Q15||Hinkel fc Kane 12013 

Duncan fc Kipping ||2Q13| |Kaltenegger||2QlQ| |Heller fc 

Armstrong 2Q14[ ^ynolds et al.||1987| ). For rocky plan- 


ets in the habitable zone, the undetected presence of a 
large satellite can confuse characterization attempts, as 
the m oon is an additi onal source of thermal or reflected 
light (Robin son 2Q11|), potent ially with its own molecu¬ 
lar features (Rem et al.j2014). Detecting large compan¬ 


ions could also aid in the characterization of potentially 
habitable exoplanets, as large satellites can help provide 
long-term obliquity, and thus climate, stability (Laskar 
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et al.||1993). However, the presence of large moons doe s 
not necessarily imply such stability (Ward et al. 2002), 
nor is a moon necessitate d for moderate obliquity stabil- 
ity (jLissauer et al.||2012|), and it might reduce habitabil- 
ity near the outer edge of the insolation habitable zone 
([Armstrong et al.||2014|). 

Given the importance of satellites to exoplanet charac¬ 
terization and in advancing our understanding of planet 
formation models, it is not surprising that many tech¬ 
niques for exomoon detection have been proposed. Most 
techniques have focused on anal yzing exoplanet transit 
timing and/or duration signals (Sartoretti fc Schneider 


19991 Szabd et al.|2006t [Simon et al.|2(JU7t|Kipping|2(JU9)' 

which can vary due to the orbit of the planet about th( 


vary due to tne orbit ot tne planet about the 

planet-moon barycenter. Additionally, mutual events, 
with a moon occulting its host star or planet while the 
planet transits the star, as well as the overlap of the 


variations in the stellar signal ([Sartoretti fc Schneider 

1999 

ISzabo et al.||2006| |Kipping||2011[ 

I'l'usnski fc Valio 

2011 

|Heller|2014[), while spectroscopy c 

luring the transit 


planet might allow for a measureme nt of the exomoon’ s 
atmospheric transmission spectrum ( Kaltenegger[[2QlQ ) 


Most of these effects have been considered in an ongoing 


se arch fo r exomoons in the Kepler dataset (Kipping et 


i 


aUplI) 


y comparison, relatively few techniques have been de¬ 
veloped that are relevant to direct observations of exo¬ 
planets, where the planet is resolved from the host star 
in either reflected or emitted light. Development of such 
techniques is prudent given the near- and long-term in- 
terest in exoplanet direct imaging and characterization 

([Kouveliotou et al. [ [2Q14| S tapel^ dt et ^ Spergel 


et al. ^15| . Gabrfca fc Schneider ( |20U7| discussed us- 

mg mutual transit and shadowing effects to detect a 
moon in an unresolved planet-moon system, which is 
made difficult by requiring an observational duty cycle 
near to 100% of the moon’s orbit and/or a fortuitous 
inclination of the moon’s orbit relative to the observer. 
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Moskovitz et al. (12009) discussed the influence of a large 
)olo 


moon on the bolomet ric therma l -infrar ed phase curve 
of its host planet, and Robinson (2011) suggested that 
phase-dependent variability at wavelengths correspond¬ 
ing to strong absorption bands for the host planet c ould 
indicate the pre sence of an airless moon. Recently, |Pe- | 


ters & Turner (2013) discussed the direct detection o: 
tidally-heated exomoons, which can be made to be very 
luminous depending on the amount of input tidal energy. 

A major obstacle to exomoon detection in direct obser¬ 
vations of exoplanets is angular resolution. For example, 
at a distance of 1.34 pc {a Cen system) the Earth-Moon 
angular separation would span 1.9 mas, and the Saturn- 
Titan angular separation at that same distance would be 
6.0 mas. Thus, to resolve these pairs of worlds, a tele¬ 
scope operating in V-band would need to have a diameter 
of about 25 m or more, and this value increases to be on 
the order of a kilometer at 10 /im in the thermal infrared. 
Here, we turn to a mature t echnique in observational 


astronomy, s pectroastrometry (Bailey 


Garcia||2C)08 ), as a new means for detecting exomoons m 


1998 Whelan & 


direct observations of exoplanets. In spectroastrometry, 
the wavelength-dependent shift in the light centroid of an 
unresolved source can be used to indicate the presence of 
an unresolved companion body and to yield orbital pa¬ 
rameters of the companion. This is limited by the optical 
design and diffraction within the telescope, which spread 
a point source of light into a broadened distribution of 
photons, referred to as the “point spread function,” or 
PSF. A spectroastrometric detection is made possible by 
the angular shift of the PSF centroid with wavelength, 
which can be measured to precisions better than the an¬ 
gular resolution of the telescope, which sets the size of 
the PSF, yielding information at the angular scale of the 
exoplanet-exomoon separation. 

We propose a new technique to detect exomoons using 
a high-contrast imaging coronagraph which suppresses 
the light of the host star exterior to an inner wor king 
angle (e.g. Gnyon et al.|2QQ6 ; Stapelfeldt et al.|2Q14 ), en¬ 


abling direct imaging of the unresolved exoplanet and ex¬ 
omoon. This can alternately be accomplished with an ex¬ 
ternal star shade ( |Cash|2Q^ ). At wavelengths where the 


exoplanet dominates the hux, the centroid will approach 
the position of the exoplanet, while at wavelengths where 
the exomoon dominates, the centroid will approach the 
position of the exomoon (Figure [T). Spectroastrometry 
then allows for the measurement m the centroid shift of 
the PSF containing the exoplanet and its moon, result¬ 
ing in the detection of the exomoon, as well as enabling 
follow-up observations for characterization of the star- 
planet-moon system. This technique works optimally for 
an instrument that is capable of measuring both spatial 
and spectral information simultaneously, although we do 
not specify the detector design in this study. 




Figure 1. Spectroastrometric signal. Left: at a bluer wavelength 
dominated by the planet, the centroid aligns with the planet’s posi¬ 
tion at the origin. Right: at a redder wavelength dominated by the 
moon (within an absorption band in which the planet’s spectrum 
is dark), the centroid shifts to the position of the moon. 

example of the characterization of a hypothetic al n earby 
planet-moon system with spectroastrometry ( ^4.2[ ). We 
discuss how our results will scale with the various as¬ 


sumed parameters (^4.3). We end with a discussion and 
conclusions. 

2. DESCRIPTION OF SPECTRAL MODELS 

We use two fiducial model systems to explore the 
prospects for spectroastrometry: an Earth-Moon twin 
around the star Alpha Gentauri A [a Gen A; 1.34 par¬ 
secs), and an Earth-Jovian analog at 1 AU around a G2V 
(sun-like) star at 10 parsecs (Table [^. We use three dis¬ 
tinct planet spectral models to simulate the spectra of 
an exo-Earth, an exo-Moon, and a warm Jupiter {i.e., 
a Jupiter-like planet at an orbital distance of 1-2 AU 
from a Sun-like star). Our Earth and Moon models are 
described in the following subsections. The spectrum 
of the warm Jupiter was ge nerated by the radia tive- 


convective models described in Burrows et al. 


Sudarsky et al. (2003). These models produce phase- 


(2004) and 


averaged spectra of irradiated gas giants for a set of 
assumed elemental abundances, and include the effects 
of cloud condensation and multiple scattering. The at¬ 
mosphere in these models is assumed to be planar, and 
the particular model that we use in this work is for a 
solar-composition, Jupiter-mass planet orbiting a G2V 
star at a distance of 1 AU. The Jovian’s atmosphere is 
taken to be in chemical equilibrium, and the equilibrium 
abundances of key atmospheric trace g ases are computed 
according to Burrows & Sharp (1999). 


2.1. Earth Spectrum Model 

To simulate Earth’s disk-integrated spectrum, we use 
the NASA Astrobiology Institute’s Virtual Planetary 
Laboratory three-dimensional, multiple scattering spec¬ 
tral Earth model, which generates temporally and spec¬ 
trally resolved disk-integrated synthetic observations of 
Earth. This model has been described and extensively 
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Table 1 

The moon and planetary parameters used in our model, where is the telescope diameter, d is the distance from the observer, tabs is 

the duration of exposure, e is the telescope efficiency factor. 


System 

Moon radius 
(m) 

Planet radius 
(m) 

Planet-moon 
separation (m) 

Orbital 
period (d) 

^tele 

(m) 

d 

(pc) 

tabs 

(hr) 

e 

Semi-major 
axis (AU) 

Earth-Moon 

1.738x10® 

6.371x10® 

3.844x10® 

27.32 

12 

1.34 

24 

0.2 

1.23 

Jovian-Earth 

6.371x10® 

6.991x10^ 

3.064x10® 

34.60 

12 

10 

24 

0.2 

1 


wavelength-dependent radiance coming from any given 
pixel is assembled from a lookup table that contains spec¬ 
tra generated over a grid of different solar and observer 
zenith and azimuth angles. Elements within the lookup 
table, which are generated using a one-dimensional, 


line-b y-line radiative transfer model (Meadows & Crisp 


1996), are computed for a variety of ditterent surface and 


atmospheric conditions, as well as several different cloud 
coverage scenarios (e.^., thick, low cloud or thin, high 
cloud). 

To simulate time-dependent changes in Earth’s spec¬ 
trum we use spatially-resolved, date-specific observations 
of key surface and atmospheric properties from Earth ob¬ 
serving satellites as input to our Earth model. Gas mix¬ 
ing ratio and/or temperat ure profiles are taken from the 
Microwave Limb Sounder (Waters et al.||2QQ6), the Tro¬ 
pospheric Emission Spectrometer ([Beer et al.|2QQl[ ), the 
Atmospheric Infrared Sounde r dAumann et al. 2UU3|) , and 
the CarbonTracker project (jEeters et al.||2U07p. Snow 
cover and sea ice data as well as cloud cover and opti¬ 
cal thickness data are taken from the Moderate Resolu- 


son et al. 

1989 

) aboard NASA's 'Terra and Aqua satel- 

lites (|Hal 

et a 

[. 1995 Riggs et al. 1999). Wavelength- 

dependeni 

t optical properties tor liquid water clouds were 


derived using a Mie theory model (Crisp 1997) and 
were parametrized usin g geometric optics loir ice clouds 
(Muinonen et al.| 1989). The spectra presented in this 
work are tor Eartn at northern vernal equinox at quadra¬ 
ture phase. 


2.2. Moon Spectrum Model 

We divide the Moon’s spectrum into two components: 
reflected solar and emitted thermal. At quadrature, 
the reflected component dominates at wavelengths below 
3.5 jim and the emitted component dominates at wave¬ 
lengths above this. The thermal component o f the spec¬ 
trum i s computed using the model presented in |Robinson| 
(2011). The model generates phase-dependent, thermal 
infrared spectra of the Moon assuming globally-averaged 
values of the lunar Bond albedo, nightside temperature, 
and the spectrally-resolved surface emissivity. 

The lunar phase function is markedly non-Lambertian, 
so we use empirical models of the Moon’s phase- 
dependent reflectivity to simulate the shortwave com¬ 
ponent of the lunar spectrum. Spectrally-resolved mea¬ 
surements of the lunar surface phase function from the 
RObotic Lunar Observatory (ROLO) are used to sim¬ 
ulate the reflected light component of the lunar spec¬ 
trum at star-planet-observer angles (z.e., phase angles) 
between 0° and 97°, which are the phase angles f or which 
the surface phase function has been published (|Buratti| 
et al. 2011). The ratio of the specific intensity emerg¬ 


ing from a surface element (/a(q^,/ io, m), where a is the 
phase angle, and /io and /i are the cosines of the inci¬ 
dent solar angle and the emission angle, respectively) to 
the incident specific s olar flux (Fy = Fq/tt) is given by 
( Chandrasekhar|p'96Q| ) , 

Ix{a,^o,p)/Fx = ^o)] , ( 1 ) 

where f{a) is the surface solar phase function (which is 
distinct from the planetary phase function), and the final 
collection of terms in brackets provides the functional 
form of the lunar scattering law. The ROLO data provide 
/(a), so that we can integrate Ix{ol^ fi) from Eq.j^over 
the illuminated disk, which eliminates the dependence on 
/io and fji^ and allows us to compute the phase-dependent, 
reflected light spectrum of the Moon, Ix{oi). 

At phase angles larger than 97°, where f{a) is not 
reported, we model the Moon’s reflect ed light spec- 
trum using the lunar phase functions of|L ane fc Irvine 
(1973), who measured the Moon’s brightness over a wide 


range of phases through a series of broadband filters 
from the near-ultraviolet to the near-infrared. By pair¬ 
ing these measurements of the Moon’s phase function 
with a medium-resolution (A/A A ^ 500) measurement of 
the disk- integrated lunar spectr um from NASA’s EPOXI 
mission (Livengood et al. |mT] i , we can infer the Moon’s 
spectruni at phase angles other than that at which the 
EPOXI data were acquired. Thus, for large phase an¬ 
gles, we take the phase-dependent, disk-integrated spe¬ 
cific brightness of the Moon to be 


Ix{a) = h{a') 




( 2 ) 


where a' is the phase angle of the EPOXI observa tions, 
and ^ is the pla netary phase function measured by [Lane 
& Irvine (1973). The lunar EPOXI observations span 
0.3 fim to about 4.5 /im in wavelength, and were taken at 
a phase angle of 75.1°. In general, our two approaches to 
simulating the reflected-light component of the Moon’s 
spectrum agree to within measurement error at phase 
angles near quadrature, a = 90°, which is the phase that 
we emphasize in this work. 


3. METHODS 

The spatial resolution needed to detect the spectroas¬ 
trometric signal for our fiducial simulated planet-moon 
systems is much greater than any telescope can currently 
provide, and greater than the proposed first-generation 
direct imaging missions. Both planet-moon systems have 
a separation of ~ 2 milliarcseconds (mas), while the Hub¬ 
ble Space Telescope has a PSE width of «100 mas at 
2.5 micron (and no coronagraph), and first-generation 
imaging telescopes will have similar spatial resolutions. 
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Consequently we focus on future second-generation tele¬ 
scopes. To simulate a spectroastrometric detection of the 
exomoon models presented in this paper, we explored the 
multidimensional parameter space of an idealized coron- 
agraphic space telescope capable of such a feat. A space 
telescope would be ideal for exomoon detection due to 
the ability to observe in the Earth’s atmospheric win¬ 
dows and to avoid atmospheric turbulence without the 
aid of adaptive optics; these same requirements apply to 
coronagraphic detection and characteriz ation of Earth- 
like planets in their habitable zones. In §3.1| we outline 
the parameters used in our mode l, in ^3.2 we define the 
spectroastrometric signal, in ^3.3 we disc uss the photon- 
limited noise and signal-to-noise, and in §3.4| we describe 
Monte Carlo simulations of spectroastrometric observa¬ 
tions. 


3.1. Model System Parameters 

We have assumed a planet-moon system whose moon 
orbits in the planet’s equatorial plane in prograde mo¬ 
tion. The planet-moon separation for the ExoEarth- 
Moon system was taken to be the current Earth-Moon 
separation, while the separation for the Jovian-Earth sys¬ 
tem was taken to be 30% of the Hill radius of a warm 
Jupiter at 1 AU from a Sun-like star (corresponding to 43 
RJupiter)^ slightly less t han the critical semi-majo r axis 
for a prograde satellite (Barnes & O’Brien 2002|). Eor 
the preliminary exploration of the signal-to-noise ratio 
for different spectral resolutions and wavelength pairs, we 
set the telescope diameter to 12 meters, distance to the 
system from the observer to 1.34 parsecs (Earth-Moon) 
and 10 pc (Jovian-Earth), and exposure time to 1 day 
(Table [^. We assume that the moon-planet position is 
fixed over the range of the exposure; since both of our 
fiducial systems have orbital periods of ^ 30 days, the 
centroid will shift by ~ 12° over this time, which is only 
20% of the planet-moon separation, and a very small 
fraction of the width of the PSE. We have selected a 
12-meter telescope based on the range of apertures be¬ 
ing cons idered for the “High-De finition Space Telescope” 


(HDST; |Dalcanton et al. 


(conserv ative) estimate o: 


2015), which is informed by a 
the detectability of Earth-like 
planets (iStark et al.||2015|). The telescope efficiency fac¬ 


tor, which accounts for instrument throughput, detector 
quantum efficiency, and photometric aperture losses, was 
kept at 20% throughout all calculations. These syst em 


parameters may be found in Table Section §4.3| ex¬ 
plores the effect that varying the telescope size cind dis¬ 
tance from the observer has on the signal to noise ratio. 

Such a large telescope would be capable of rapid detec¬ 
tion of Earth-sized exoplanets (in of order a few weeks) 
with a dire ct-imaging survey of nearby Sun-like stars 
(|Agol||2007|). The planets found with this initial coron- 
agraphic imaging survey could then be prioritized based 
upon, for example, their presence within their host-star 
habitable zone and potential for the existence of a large, 
detectable exomoon, to motivate the additional observ¬ 
ing time for detailed spectral and spectroastrometric 
characterization. 

We note that the parameters of the Jovian-Earth 
model system we have chosen are driven by favorable 
observational detectability rather than theoretical preju¬ 
dice. Models of satellite formation around giant planets 
have been tailored for the Solar System to reproduce, for 


example, the mass ratio of 1 0 ^ observed for the satellite 
systems of t he giant planets (jCanup & Ward|2002| | Sasaki 
et^lMol ), SO tha t giant planets may harbo r satellites 


as large as 0.7 (Heller & Pudritz 2015b). In addi¬ 
tion, models of in-situ formation of satellites from a disk 
produce regular satellites at distances of 6 — SORjupUer- 
Thus, within the context of these models, an Earth-sized 
satellite at 43 Rjupiter distance from a giant planet is 
unexpected, but may have formed by other means, in¬ 
cluding capture and migration. 


3.2. Speetroastrometrie Signal 

The spectroastrometric signal is the difference in the 
centroid positions when the system is observed at two 
different wavelengths. We define the centroid position, 
c(A), as being the fiux-weighted center of light of the 
planet and moon as a function of wavelength. A, in the 
direction of two angular sky coordinates, (x, y), which we 
denote as c(A) = (ca^(A), Cy(A)), where Cx,y are the {x^y) 
centroids with respect to the origin of the star-planet- 
moon system (and the orientation of the coordinates is 
chosen with respect to an observer-defined reference di¬ 
rection). Note that the angular coordinates (x^y) are 
expressed in milliarcseconds (mas), which is the typical 
order of the moon-planet separation. 

The detector will be capable of measuring this signal 
over a wave band with central wavelength A and spectral 
resolution R, giving the width of the band A A = X/R. 
We define the band-averaged centroid as 


c(A) 




r X-\- j 2 
Ja-AA/2 


XFxdX 


( 3 ) 


where F\ is the total fiux of the planet and moon, and 
we have weighted the centroid by the number of photons 
detected. 

In computing the model centroid, we define and 
lo represent the fiux of the moon and the planet, 
respectively, and and Vp to represent the positions of 
these bodies (in two dimensions, projected onto the sky 
plane, perpendicular to our line of sight). The position 
of the wavelength-dependent centroid, c(A), is 


Rm,X^m F R P,AI*p 
d{Fm,x F Rp,x) 
Rm,xP ^p 

Rm,X F Fp^x ^ 


( 4 ) 


where p = (r^ — Yp)/d\^ the angular sky separation be¬ 
tween the moon and planet in radians at a distance d 
from the observer (we convert this to milliarcseoconds). 
Note that this equation neglects the variation of the cen¬ 
ter of light of each body with the illumination phase; this 
is on the order of the angular size of each body, which is 
much smaller than the spectroastrometric signal. 

Measuring c(A) requires a reference position on the sky, 
for example the planet’s position, while the planet’s po¬ 
sition is not known a priori in the presence of a centroid- 
shifting moon. However, knowledge of the planet’s po¬ 
sition is not needed when considering only the differ¬ 
ence between the centroids measured in different spectral 
bands. This is the spectroastrometric signal, 5'(Ai,A2), 
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which scales as the scalar difference of the planet/moon 
centroid in two wave bands, 

5(Ai,A2) = |c(Ai)-c(A2)|, (5) 

which we measure in milliarcseconds. Note that it may 
be possible to select a different resolution R for each 
wave band. A significant advantage of measuring the 
spectroastrometric signal, 5'(Ai,A2), rather than c(A) is 
that the absolute position of the planet-moon system does 
not need to be calibrated; only the relative position with 
wavelength needs to be measured. This obviates the need 
for an absolute sky coordinate reference frame, which 
can be a challenging measurement to make (although 
the host star could perhaps be used, as discussed below). 
Also, the position of the planet is affected by its orbit 
about the star, its illumination, its orbit about the center 
of mass with the moon, and its gravitational perturbation 
by companion planets; measurement of the spectroastro¬ 
metric signal eliminates all of these confounding effects. 
Note that this technique does not require flux calibration 
since the spectroastrometric signal is normalized by the 
total flux of the planet and moon. 


3.3. Signal to Noise Ratio of Deteetion 

The spectroastrometric noise scales as the ratio of the 
width of the point spread function incident on the detec¬ 
tor to the square root of the number of photons, assum¬ 
ing Poisson noise from the planet and moon dominates 
the uncertainty. In light of the extraordinary engineering 
called for in the success of this project, we have assumed 
that an ideal coronagraphic space telescope would pro¬ 
duce negligible instrumental noise (due to the dark cur¬ 
rent, read noise, scattered light, pixel size, etc). We esti¬ 
mate the number of photons, N{X), incident on the coro- 
nagraph’s detector due to the planet and moon within a 
band of width AA for some exposure time tobs to be 


r)2 ^A-|-AA/2 \ 

N{X) = iT€^tobs / j-FxdX 

4 Jx-AX /2 he 

AA A^ 


^ 7re 


r)2 

-^tele ‘ 


A he 


l^obs^X-) 


(6) 


where e is the chosen efficiency factor of the telescope 
(defined as the fraction of all photons incident on the 
telescope aperture that are measured by the detector; 
we assume that it is wavelength-independent), is 

the diameter of the telescope, h the Planck constant, c 
the speed of light, and Fx is the wavelength dependent 
flux density (the approximation assumes that this flux 

is roughly constant across a narrow bandb _ A deepe r 

analysis of this calculation can be found in Agol (2007); 
note that R = A/AA is the spectral resolving power. 

We assume that the PSF of the planet is well approxi¬ 
mated by an Airy disk, which we in turn approximate as a 
Gaussian with angular profile I{0) oc eyeg/(T psf)‘^^ 


where (JpsF = 0.45A/DteZe is the angular size of the Airy 
disk, and 0 is the angular coordinate from the center of 
the PSF. The precision of the measurement of the cen¬ 
troid improves with the square root of the number of 
photons detected. Thus, the spectroastrometric noise for 
an ideal (photon-noise limited) coronagraphic space tele¬ 


scope is then 

a PSF 


0.45A 


A(A)V2 D 


7re 


4 A ^ 


cr(A) = 

UtpJ.e \ ^ A fiC / 

( 7 ) 

The spectroastrometric signal-to-noise ratio (hereafter 
SNR) we define as 


SNR = 


g(Ai,A2) 

^/FXlF+~a{XF^ 


(8) 


between two wavelengths Ai and A 2 . Below we set the 
threshold for exomoon detection in our model systems at 
a minimum 5—cr confidence level, i.e. SNR > 5. 

In practice, the noise is typically dominated by the 
moon wave band. This gives an approximate expression 
for the signal-to-noise ratio of 


SNR^ 


h^tele 

O. 45 A 1 



AAi A? 
Ai he 


tobs^^Xi 


1/2 


( 9 ) 


where we have assumed that Ai is dominated by the 
moon while A 2 is dominated by the planet, and ^ = \/3\ 
is the moon-planet separation on the sky at the time of 
observation. Note that in the wave band dominated by 
the moon, Fx oc (R^/d)^, so overall the signal-to-noise 
scales as 


SNR oc ajYipRjYi 


h^tele 

d 


R ) \nRl, 


1/2 


( 10 ) 

where Omp = |r^ — r^j = /3d is the sky-plane projected 
physical separation of the planet and moon, R^ is the 
radius of the moon, and the last quantity in this equation 
is the disk-integrated specific intensity of the moon. 

We use these equations below for estimation of the de¬ 
tectability of a moon-planet system, while for charac¬ 
terization of a detected system, we employ Monte Carlo 
simulations, described next. 


3.4. Monte Carlo Simulations 

We have carried out Monte Carlo simulations of coro¬ 
nagraphic observations of a detected moon-planet sys¬ 
tem with the same idealized assumptions: we assume 
a broad range of wavelength sensitivity, and we assume 
that the photon counting noise and diffraction limit can 
be achieved. We define spectral bins equally spaced in 
log wavelength between a minimum and maximum wave¬ 
length. For each spectral bin, we compute the predicted 
number of photons within the bin for each body, A^(A), 
and then draw the observed number of photons from a 
Poisson distribution. We next add a random Gaussian 
deviate in both directions on the sky with a standard de¬ 
viation of cr(A) = 0.45A/(DteZeA^(A)^/^) at the position 
of each body. We then compute the flux-weighted aver¬ 
age positions of both bodies to obtain c(A), and compute 
the centroid uncertainty, cr(A), from a weighted mean of 
the individual bodies’ centroid uncertainties. The result 
of this is a total photon flux and (x, y) centroid position 
for each spectral bin, c(A), as well as the photon-noise 
limited uncertainty (which we assume to be the same in 
the X— and y— directions). 

We have also run simulations in which we randomly 
draw individual photons from the spectral shapes of each 
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body, and assign a position and positional uncertainty to 
each photon. We then bin these photons in wavelength to 
measure the spectrum, and compute the centroid within 
each bin to measure the spectroastrometric signal, c(A); 
we compute uncertainties on the photon flux from the 
square root of the number of photons in each bin, and 
on the centroid from the scatter in the photon positions 
divided by the square root of the number of photons in 
each bin. This approach is more time intensive, yet gives 
equivalent results to the pre-binned spectrum approach. 

4. RESULTS 

Our results outline the range of instrumental, spec¬ 
troscopic, and planetary system parameters for which 
we expect the spectroastro metr ic method to detect the 
presence of an exomoon. In §4.1[ we discuss the optimum 
observing wavelengths and spectral resolutions found 
with telescope and planet- moon syst em parameters 
fixed at values described in §3.1| In §4.2| we give an 
example of characterizatio n of an Earth-Moon twin at 
the distance of a Cen. In §4.3[ we discuss the results of 
varying the distance and telescope size for both systems 
while spectral resolution and observing wavelengths are 
fixed at favorable values. 


4.1. Optimum Spectral Resolutions and Observing 
Wavelengths for the Model Systems 

Sections §4.1.1| and §4.1.2 discuss our particular model 
systems in more detam Results for both systems re¬ 
flect the inherent challenge of spectroastrometry: moon- 
dominated wavelengths tend to exist because the host 
planet is dim due to atmospheric absorption, not due to 
emission features of the moon. Thus, the bands with 
highest centroid offset inherently have less total flux and 
an elevated level of noise. In order to bring the SNR at 
these wavelengths above the detection threshold, a bal¬ 
ance must be struck between a spectral resolution which 
is low enough to let in many photons, yet high enough 
to reveal a measurable shift in the centroid. 


4.1.1. Moon-like Exomoon Orbiting an Earth-like Planet 

In an Earth-Moon analog system, the Moon outshines 
the Earth in the water bands around 1.96 and 2.6 - 3.0 
/rm, the carbon dioxide band at 4.2 - 4.4 /rm, and it 
has a sizable thermal excess around 5.0 - 8.0 jam (Eig. 

top plot), although in this paper we only consider 
wavelengths out to 3 micron. The maximum fraction 
of the Moon’s flux for the Earth-Moon system occurs at 
A = 2.69 /im, contributing 99.8% to the total flux when 
viewed with spectral resolving power R = 100 (Eig. 
bottom plot). 

To optimize the spectroastrometric SNR requires 
choosing two wavelengths and resolutions at which the 
difference in centroid is large, but the noise is also small. 
At lower spectral resolution, the centroid offset, c(A), be¬ 
comes smoother so that the molecular absorption bands, 
at which the Moon dominates the flux, are avera^d with 
wavelengths where the Earth dominates (Eig. 1^. The 
noise decreases at lower spectral resolution due to the 
larger number of photons in each band, which compen¬ 
sates for the smoothing of the centroid offset. Conse¬ 
quently lower spectral resolution generally gives a higher 
signal-to-noise of detection. 


Our calculated SNR results for the Earth-Moon system 
at 1.34 pc exceed the 5—a detection threshold for the pa¬ 
rameters listed in Tabled for a wide range of wavelengths 
and resolutions (Tabled Eig. [^. The optimal wave¬ 
lengths cluster near 0.35 micron and 2.7 micron where 
the planet and the moon dominate the flux, respectively. 

4.1.2. Earth-like Exomoon Orbiting a Jovian 

In a Jovian-Earth system, the Earth-like moon out¬ 
shines the Jovian in the NIR methane absorption bands, 
shown in Eigurej^ The maximum fraction of the moon’s 
flux for the Jovian-Earth system occurs at A = 1.83 /rm, 
contributing 99.1% to the total flux, although the use of 
this band in the spectroastrometric signal does not pro¬ 
duce the highest SNR due to its higher noise; a higher 
SNR is achieved in the methane band near 1.4 /rm, albeit 
at slightly smaller centroid offset. 

The Jovian-Earth system at 10 pc produces maximum 
SNRs above 5—a for the parameters in Table assum¬ 
ing spectral resolutions between R = 2 to R = 500, and 
wavelength pairs drawn from 0.430 - 0.571 fim and 0.858 
- 0.947 jam. Table lists the optimum wavelength pairs 
for a selection of spectral resolutions. To illustrate how 
the SNR changes as a function of spectral resolution, 
Eigurej^ displays the SNR for a range of combinations of 
wavelengths for two different spectral resolutions. Note 
that in this model system several wavelength pairs would 
achieve a SNR well above the detection threshold of 
5—cr, although the optimum resolution will ultimately 
depend upon the instrument design. Eigure shows a 
slice through this plot with the choice of one wave band 
centered at 0.89 /rm; this shows that a broad range of 
comparison wavelengths could be used. 


4.2. Characterization of an Earth-Moon system around 

a Cen A 

Once a planet-moon system is identified with spec¬ 
troastrometry, more telescope time will be expended on 
a habitable-zone system with astrobiological interest in 
an attempt to characterize the moon-planet system more 
fully. If spectroastrometric measurements can be made 
simultaneous to spectroscopic measurements, then the 
time will concurrently be used for high precision spec¬ 
troscopy to characterize the atmosphere and look for 
molecular signatures, as well as to more precisely char¬ 
acterize the spectroastrometric signal. 

In particular, the spectroastrometric offset can be used 
to separate the individual spectra of the planet and 
moon, to measure the orbit of the moon, and to mea¬ 
sure the mass of the planet. As a concrete example, we 
consider an Earth-Moon twin system observed for a total 
exposure time of one month. The motion of the moon 
about the planet will cause the centroid to vary with 
time. In the case of a circular, face-on system, the angu¬ 
lar motion of the centroid can be followed with time to 
trace out the (nearly) Keplerian orbit of the moon about 
the planet. Once an orbit is completed, the angular mo¬ 
tion can then be removed to give the centroid offset as a 
function of time; then, as a function of wavelength this 
can be binned to give the scalar centr oid o ffset. 

We ran a Monte Carlo simulation (^3.4) of an Earth- 
Moon twin orbiting a Cen A, which is the most favor¬ 
able case for detection and characterization of a moon 
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Figure 2. The flux of a twin Earth-Moon system at quadrature phase angle orbiting a Cen A as a function of wavelength in microns 
(top) and the contribution of flux due to the moon shown as a fraction of the total flux (bottom). The maximum fraction of the moon’s 
flux for the Earth-Moon system occurs at A = 2.7 /j,m, contributing 99.8% to the total flux. 


with spectroastrometry; the proximity of t his star fa¬ 
vors many modes of exoplanet detection ( |Eggl et ah 


2013). We assume that the planet/moon instellation is 


the same as Earth’s, which puts the semi-major axis at 
1.23 AU; planet orbits can be stable for l ong timescales at 


this distance given the right properties (Andrade-Ines & 


Michtchenko 12014 ). We ignore the effect of the compan¬ 


ion star {a Cen B) on the observation; however, this may 
contribute a potential source of noise. Since the effective 
temperature of the star {T^ff = 5790 K) is slightly hot¬ 
ter than the Sun, we multiply our simulated spectra by 
the ratio of Planck functions at the two temperatures, 
divided by the ratio of the temperatures to the fourth 
power (to maintain the same incident flux). 

In Eigure we used the Monte Carlo simulations of 
a 28-day direct imaging observation of the Earth-Moon 
analog system (orbiting a Cen A) to measure the cen¬ 
troid offset, assuming that the motion is face-on, circular, 
and that the angular motion can be corrected for. This 
shows that a SNR can be achieved that varies with wave¬ 
length between the extremes of the planet centroid (at 
zero) and the moon centroid (at 1.9 mas). The detec¬ 
tion is highly significant, a total of 169—a (in practice 
additional sources of noise may reduce this value). 


4 . 2 . 1 . Spectral Distentanglement 


We can then use this measurement to approximately 
separate the spectra of the planet and moon as follows. If 
we assume that the extremes of the centroid motion with 
wavelength correspond with the positions of the planet 
and the moon, then we can identify the two wavelengths 
at these extremes, Ai and A 2 respectively, and compute 
the centroid offset signal, Smax = 5'(Ai,A2). We can 
then divide the centroid offset relative to the planet- 
dominant wavelength, 5'(A,A2), by the maximum offset, 
Smax- Then, since the center of light varies between the 
planet and the moon, the fraction of flux from the moon 
is simply 

fmoon — ^ 2 )/Smax-i (H) 

and the planet is 

fplanet = 1 “ S{X, X2)/Smax- (12) 

We can then multiply this by the total flux, Ffoti fo 
recover the spectra of the individual bodies, Fmoon = 
fmoonFtot and Fplanet = fpianetFtot- An example of this 
based on the simulated 28-day observation of the Earth- 
Moon twin orbiting a Cen A (at d = 1.34 pc) with a 
12-meter telescope and 300 wavelength bins from 0.3 to 
3 micron {R = 130) is shown in Eigure The recovered 
spectra of the bodies match the input spectra, and this 
approach allows ‘cleaning’ the spectrum of the planet 
from the contribution from the moon. This could enable 
spectral retrieval to be carried out on the individual bod- 
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Figure 3. The centroid offset, c(A), and uncertainty, cr(A) for an Earth-Moon twin at the distance of 1.34 pc (aCen; see Table [n for 
simulation parameters). The reference position is the centroid of the planet (c(A) = 0 mas), while the moon is at c(A) = 1.9 mas (red 
dotted line) and the center of mass is at 0.023 mas (green dotted line). Black lines are c(A), while blue are cr(A). Solid, dotted, and dashed 
black/blue lines correspond to R = 100,10,1. 


ies, constraining their atmospheric properties and com¬ 
positions with further modeling. 

If it cannot be determined whether the moon-centered 
wavelength is truly dominated by the moon, then an¬ 
other approach is to estimate the semi-major axis from 
the mass measured with another technique, such as re¬ 
flex motion of the star measured with astrometry or ra¬ 
dial velocity. Then the predicted semi-major axis can 
be computed from the period of the spectroastrometric 
variation, and then this may be used to disentangle the 
spectra assuming one wavelength is dominated by the 
planet. 

4.2.2. Planet Mass and Orbit Characterization 

The time dependence of the spectroastrometric signal 
can be used to measure the mass of the planet. Based 
on the simulation of a 28-day signal for a face-on Earth- 
Moon system at 1.34 pc, we chose the wavelength bins 
dominated by the Earth and Moon, respectively at 0.35 
and 2.75 micron with R = 130. We then measured the 
difference between the centroids at these wavelengths 
with time assuming 1-day exposures; Eigure shows 
a comparison of the input positions, the measured posi¬ 
tions, and the best-fit circular orbit (this can be easily ex¬ 
tended to an elliptical or edge-on orbit, but we picked this 
face-on system as an example). We ran 100 realizations 
of the observations, and we found M = 1.03 ± O.I 2 M 0 
(the measured mass includes the mass of the moon). This 
relies on having two wavelengths at which the moon and 
the planet dominate the flux; it may be possible to im¬ 
prove the precision with modeling of the data, which we 


leave to future investigation. 

4.3. Distance and Telescope Size Dependence 

The SNR’s simultaneous dependence on distance and 
telescope diameter for the Jovian-Earth system is shown 
in Eigure IE Using the wavelength pair Ai = 0.45 /am 
and A 2 = 0.89 /am at R = 50 and tobs = 24 hrs, a 6.7- 
m telescope would be sufficient to detect an Earth-like 
moon orbiting a Jovian around a Sun-like star at 10 pc 
using spectroastrometry with a confidence of 5—cr. Ac¬ 
cording to the Gliese Catalog of Nearby Stars, there 
are about 20 main sequence E, G, and K type stars 
within 10 parsecs that could be used as potential tar¬ 
gets. These systems and the diameters needed to detect 
such a planet-moon system with a confidence of 5—cr may 
be seen in Eigure 12 Since the effective temperature of 
each star differs from the Sun, we multiply our simulated 
spectra by the ratio of Planck functions at the two tem¬ 
peratures, divided by the ratio of the temperatures to the 
fourth power (to maintain the same incident flux, which 
is assumed to be the Solar constant). This assumes that 
the albedos of the planets remain the same for planets 
orbiting stars of different spectral types, and that the 
stellar spectral features differ slightly; both assumptions 
are approximate, but we expect will have a small effect 
on our results compared to, for example, planets with 
different atmospheric compositions that what we have 
simulated. 

The SNR’s dependence on distance from the observer 
arises from both the signal’s inverse relationship to the 
distance and the flux’s contribution to the noise, where 
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Figure 4. The spectroastrometric SNR of the Earth-Moon system shown for a combination of wavelengths at R = 1.5 (top) and R = 5 
(bottom). Any point on the plots corresponds to a pair of wavelength bands marked on the x and y axes. The color contours indicate the 
level of SNR achieved when the centroids of those two wavelength bands are compared. We assume a SNR of 5—cr or greater would yield 
a detection. All calculations were done for a model system analogous to our Earth and Moon located 1.34 parsecs away in the habitable 
zone of a Cen A, with an exposure time lasting 24 hours, and a 12-m telescope diameter. 


the total flux Fx is inversely proportional to the square of 
the distance d. Equation 10 shows that the SNR scales as 
where one power comes from the decreasing angular 
separation of the moon and planet, and one power comes 
from the noise caused by decreasing flux of the moon. 

The telescope diameter affects the width of the point 
spread function as as well as the total number 

of photons incident on the detector in a given exposure 
time fobs as giving a scaling of the SNR as as 


shown in equation pR| 
SNR as (Aeie/rf)^ tK 


Given the overall dependence of 
e volume which can be surveyed to 


a fixed SNR scales as oc (holding other parame¬ 


ters fixed). If the diameter of the telescope is decreased, 
then the observing time needs to increase as fobs oc 
to obtain an equivalent SNR; thus, either fewer systems 
can be surveyed due to the longer observing time, or the 
observing time will increase until systematic errors dom¬ 
inate the noise. 

We have focused on habitable-zone planet/moon sys¬ 
tems around G stars in this paper, but moons should 
be searched for around planets at a range of star- 
planet/moon separations. More distant planet-moon 
pairs will have lower intensity due to the decreased flux 
they receive from the star: Fx oc a“^, where a is the semi- 














10 


Agol, Jansen, Lacy, Robinson & Meadows 




Figure 5. The flux of the Earth-like moon and the warm Jupiter at quadrature phase angle as a function of wavelength in microns (top) 
and the contribution of flux due to the moon shown as a fraction of the total flux (bottom). The maximum fraction of the moon’s flux for 
the Jovian-Earth system occurs at A = 1.83 /rm, contributing 99.1% to the total flux. 


major axis of the orbit of the planet and moon about the 
star. However, their Hill spheres will expand due to their 
larger separation from the star, so if the moon’s orbital 
distance is a fixed fraction of the Hill sphere, then the 
spectroastrometric signal will increase oc a due to the 
wider moon-planet separation. These two effects com¬ 
pensate to cause an equal SNR as a function of the dis¬ 
tance of the planet/moon system from the host star. 

5. DISCUSSION 


The detection of an exomoon is a long-sought discov¬ 
ery which could profoundly impact our understanding of 
planet and satellite formation and evolution, including 
that of our own Solar System. It is clear that moons 
can outshine their planets at some waveleng ths and con¬ 
tribute to apparent variability of a planet ([Williams fc 


triDute to apparent variaDility oi a planet »W iJ 

Knacke|[2QQ4 Moskovitz et al. 2009 [Robinsrap 


m 

practice, though, it may be difiicult to prove that spec- 
troscopic or time-dependent signatures are in fact due to 
a moon. Since a moon should be spatially separated from 
a planet, and the separation should follow a (nearly) Kep- 
lerian orbit as a function of time, spectroastrometry gives 
a means of detection that is more definitive, and would 
allow measurement of the properties of the planet-moon 


system. 

If the separation of the planet and moon is similar 
to the width of the instrumental point-spread function 
(PSF), then one can simply resolve the planet and moon. 
For example, a favorable configuration of a satellite or¬ 
biting at 0.25 of the Hill radius of a Jovian-mass planet 
orbiting at ^ 1 AU around a star at 3 pc will have a 
maximum angular separation of ~ 6 mas. This is about 
the angular resolution at 0.5 /im of an 8-meter space tele¬ 
scope; in practice, the angular separation will be smaller, 
so most planet-moon systems will be unresolved, hence 
the need for spectroastrometry. The a Cen A case of an 
Earth-Moon analog can almost be resolved directly: at 
0.35 /im, a 12-meter telescope has a resolution of ~ 2.7 
mas, while the angular separation of the Earth-Moon 
analog is 1.9 mas. Since the Moon is much fainter than 
the Earth at these wavelengths (< 1%), it may be diffi¬ 
cult to deconvolve the light, and thus spectroastrometry 
will still be favored for detection of the moon. 

We have chosen 3 micron as the long-wavelength limit 
for these observations based on the fact that current coro- 
nagraphic designs are considering 2 micron as a possible 
wavelength cutoff, with a possible extension to 3 - 5 mi- 
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Figure 6. The spectroastrometric SNR of the Jovian-Earth system shown for a combination of wavelengths at R = 4 (top) and R = 50 
(bottom). Any point on the plots corresponds to a pair of wavelength bands marked on the x and y axes. The color contours indicate the 
level of SNR achieved when the centroids of those two wavelength bands are compared. We assume that a SNR at or above 5—cr could 
yield a detection. All calculations were done for a model Jovian-Earth system located 10 parsecs away, orbiting at 1 AU around a G2V 
star, with an exposure time lasting 24 hours, and a 12-m telescope diameter. 


cron (Dalcanton et al. 2015); longer wavelengths have 
the problem of a larger inner working angle and noise 
due to thermal emission from the telescope. Spectral 
coverage that extends to 3 /im covers the water band at 
2.7 /im and the methane band at 2.3 /im; these bands 
have the advantage of being dominated by the moons 
in the two cases we have considered. In particular, for 
measuring the mass of the Earth-like planet in the Earth- 
Moon analog case, the 2.7 /im band is critical for being 
able to measure the planet-moon semi-major axis. At 


3 /im we find that there are ~ 7 KV stars out to 5 par¬ 
secs which can be resolved at > 2A/DteZe for = 12 

meters; 15 GV stars out to 17 parsecs; and 36 EV stars 
out to 24 parsecs (we have not considered stars beyond 25 
parsecs). The possibility of detection of an exomoon in 
these systems depends strongly on the potential proper¬ 
ties of the exomoon-exoplanet systems present, including 
their frequency. 

Another issue we have not considered is the effect on 
these observations of exo-zodiacal light due to scatter- 
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Figure 7. Methane absorption in the Jovian atmosphere leads to the moon-dominated band with a central wavelength of 0.891 [im for 
R = 50 (marked for reference on figure by blue dashed vertical line and shaded green region). Fixing the centroid from this band as one 
of our comparisons, the plot shows the SNR (plotted as a black line) for each other band with resolving power R = 50 where signal is 
measured relative to that fixed 0.891 iim band. 



1 

Wavelength (micron) 

Figure 8. (Top) Centroid offset, c(A), versus wavelength for a simulated observation of an Earth-Moon twin at 1.34 pc. The orbital 
motion of the system has been removed (it is assumed to be face-on). The black line is along the axis connecting the Earth and Moon; 
the red line is perpendicular. The zero point is centered on the Earth. (Bottom) Significance of the centroid offset relative to the mean 
centroid in the directions parallel (black) and perpendicular (red) to the line connecting the Earth and Moon. 
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Figure 9. Recovered spectra of the Earth (black) and Moon (red) using the spectroastrometric offset. The lines show the input spectra; 
the dots show the recovered spectra. 
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Figure 10. Simulated spectroastrometric orbit of the Moon-Earth system at 1.34 pc, which is the difference between the A = 2.76 and 
0.35 micron with R = 130. The red curve shows the input positions, the crosses show the measured positions at 1-day intervals, while the 
blue shows the best-fit ellipse for one Monte Carlo realization of the observations. 
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Figure 11. How the SNR scales with telescope diameter and distance from the observer for the Jovian-Earth system observed between 
0.45 /rm and 0.89 /rm at R = 50. The white region shows SNR > 70—cr, and the grayscale shows varying SNR strengths, where the black 
region shows a non-detection (<5—cr). The green line indicates the telescope diameter needed to directly resolve the Earth-like exomoon 
from the Jovian planet at 0.89 fim as a function of distance, in which case spectroastrometry would not be the primary method to detect 
an exomoon. 


ing/thermal emission by dust. This diffuse emission com¬ 
ponent could have a spectrum similar to that of a moon, 
and if concentrated in a region near the planet, may be 
confused with the effect of a moon. If spread throughout 
the system it would contribute to the noise of the obser¬ 
vation, which grows more severe at longer wavelengths 
due to the larger PSF, which scales in solid angle (or 
area on the detector) as A^. If (exo-) zodiacal light is a 
significant source of noise , then shor ter wavelengths may 
be favored. We estimate ( Agol|2007|) the ratio of the flux 
of the exozodi flux within the planet’s PSF, Fez^ to flux 
of the host star, to be: 


Fez cP „ / a 

-^ = ^TEzSf 


o ' H/Z/J CLC I / 

^ \Ftele J 


-2 

(13) 


where tez (units sr“^) relates the incident stellar flux to 
the exozodi surface brightness, and r is the planet-star 
separation. This is less than the Moon flux for Solar- 
System like exo-zodiacal light, but it can increase with 
distance, wavelength, or dust surface density, and thus 
may play a factor for more distant, young systems. Older 
systems may be largely void of dust (like our Solar Sys¬ 
tem). 

The precise choice of the wave bands used for detection 
will depend on multiple factors rather than just optimiz¬ 
ing the signal-to-noise ratio. For example, if a telescope 
has a short or long-wavelength cutoff that differs from 


what we assumed, this may require a different choice of 
band. If a narrow wavelength range is required, then it 
may be necessary to choose more closely spaced wave¬ 
lengths than what we have chosen. As an example, con¬ 
sider the fiducial Earth-Moon twin system with parame¬ 
ters from Table We tried all pairs of wavelengths from 
0.3 - 3 micron, and resolutions from R=I-I0 (in steps of 
I) which give SNR > 10 for this simulated system (for a 
12-meter telescope and I-day observation). From these, 
we found wave bands that give the largest spectroastro- 
metric signal, 5'(Ai, A 2 ): we find that Ai, A 2 = 0.335, 2.73 
micron and R= 10, gives a centroid offset that is 90% of 
the Moon-Earth sky separation at quadrature (SNR = 
14). Higher resolutions can approach 100% (Eig. 1^, but 
at lower SNR; this is wha t wa s required for the Sarac- 
terization of the system (^4.2). 

Another consideration is the separation in wavelength: 
if the instrumental systematic errors of the centroid mea¬ 
surements grow with the separation in wavelength, then 
it may be advantageous to choose two bands that are 
closely spaced in wavelength. Amongst the SNR >10 
band pairs for the Earth-Moon system at 1.34 pc, we 
find Ai,A 2 = 2.3,2.7 micron (R=I0) gives a centroid 
offset that is 75% of the sky separation. The precise 
choice of bands and resolutions will depend strongly on 
the planet properties and stellar distance, as well as the 
properties of the coronagraphic telescope, so these exam¬ 
ples are merely for illustration. 

With the detection of a moon with spectroastrometry, 
there are several applications of the detection which may 
be used for characterizing the planet-moon system. A 
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spectroastrometric detection will confirm that a spectro¬ 
scopic moon signal is truly due to an exomoon, and not, 
for example, due to varying temperature or gas mixing 
ratio profiles within the planet’s atmosphere; it will also 
confirm that the planet is in fact a planet. The variation 
of the spectroastrometric signal with time will follow the 
orbit of the moon about the planet, giving the period of 
the moon’s orbit when monitored with time. As shown 


above (^4.2), using Kepler’s law, the period and semi¬ 
major axis allow the determination of the mass of the 
planet plus moon; this could enable the mass measure¬ 
ment of an Earth-Moon twin at 1.34 pc with a 12-meter 
telescope at 0.35 and 2.76 micron. This will also po¬ 
tentially allow the determination of the semi-major axis, 
orbital inclination, and eccentricity of the moon, if the 
signal-to-noise is sufficient. These applications require 
that wavelengths can be identified at which the planet 
and moon each dominate the flux, although this may be 
difficult to prove in practice as the spectroastrometric 
signal is degenerate between the moon’s radius, and 
semi-major axis, S oc R‘^amp- 
In addition to the spectroastrometric signal, which 
yields the difference of the center-of-light between the 
moon and the planet, if the absolute astrometric signal 
of the planet can be measured, then it might be possible 
to also measure the mass of the moon due to the reflex 
motion of the planet about the center of mass with its 
moon. However, this signal, which is of order the size of 
the Earth for the Earth-Moon system, may be swamped 
by the contribution of the Moon’s signal to the centroid, 
which would need to be removed precisely to measure 
the offset of the Earth from the center of mass with the 
Moon; this will require identification of the wavelength at 
which the Moon-Earth flux ratio is less than the Moon- 
Earth mass ratio, such as near 0.35 micron (Eig.J^. In 
addition, the center-of-mass offset of the planet’s light is 
in the opposite direction of the of the cent er-of-light off¬ 
set (relative to the center of mass) at wavelengths where 
the moon has a larger flux fraction. There is a wave¬ 
length at which the two offsets cancel, at which point the 
center of light coincides with the center-of-mass; for the 
Earth-Moon system at quadrature this occurs at 0.493 
/im (Eigurej^. The astrometric signal will also be af¬ 
fected by the non-uniform illumination and non-uniform 
albedo (or thermal emission) of the Earth, which will 
vary with orbital phase and rotation of the planet, which 
should be accounted for if the astrometric precision is 
sufficiently high. This could, in fact, provide another 
opportunity to detect spatial variations in the surface of 
the planet, and potentially reveal the full obliquity of the 
planet, recirculation of heat (which would be affected by 
the atmospheric pressure and opacity of the planet), and 
surface markings (or persistent cloud features). In ad¬ 
dition, perturbations by other planets in the system can 
affect the motion of the center-of-mass; for the Earth- 
Moon system, these perturbations also have an ampli¬ 
tude that is comparable to the radius of the Earth. De¬ 
spite these complications, this measurement could also 
be attempted with a coronagraphic telescope capable of 
precision astrometry of planets, which could be measured 
with respect to the position of the host star. A diffractive 
element introduced into a high-contrast imaging corona- 
graph can yield the astrometric position of the star, and 
thus allow the precise measurement of the planet’s orbit 


with respect to the star (Guyon et al.||20I3). 

A spectroastrometric detection of an exomoon can po¬ 
tentially result in a detection of eclipses/transits/mutual 
events. The orbit of the moon about the planet mea¬ 
sured from spectroastrometry may allow one to forecast 
whether and when a moon/planet will transit. Knowing 
whether the events will occur and observing at the fore¬ 
cast times gives a significant advantage over trying to ini¬ 
tially detect an exomoon with eclipses/transits/mutual 
events for which continuous observations are required 
and a successful dete ction depends on the unkn own ge¬ 
ometry of the orbits (Cabrera & Schneider||2007 ). There 
are four possibilities o f the system geometry to consider 
(Schneider et al. pM5| : L) if the moon’s orbit has a small 
inclination with respect to the planet’s orbit (which can 
be measured from its motion about the star with time), 
then the moon can regularly pass between the star and 
the planet, causing a shadow to be cast on the planet, 
which may be visible by a dimming of the planet at the 
predicted time of the stellar eclipse; 2) the passage of 
the moon into the planetary shadow (with respect to the 
star) can also yield a dimming due to the lunar eclipse; 3) 
if the inclination of the moon relative to our line of sight 
is nearly edge-on, then one might observe a transit of the 
plane t by the moon, causing a dimming at the predicted 
time (jLivengood et al.||20II ); 4) likewise, the moon can 
pass behind the planet. It is possible for cases 1/2 and 
3/4 to occur simultaneously (given a fortuitous geome¬ 
try), causing an even larger depth. The latter two cases 
have a ~ 2% probability for a randomly placed observer 
of the Earth-Moon system, which is about four times 
larger than the probability of viewing the Earth to tran¬ 
sit the Sun (^ 0.5%). Any four of these configurations, if 
detected at sufficient signal-to-noise, might yield a mea¬ 
surement of the radii of the planet and/or the moon (or 
their ratio), giving a means of determining the density 
of the planet (and possibly the moon). In addition, if 
the spectral dependence of the transits/eclipses are mea¬ 
sured, this might yield th e detection of spatially-res olved 
features on either body. Heller & Albrecht | ([2014) have 
also studied the possibility of measuring the Kossiter- 
McLaughlin effect in high-resolution spectroscopy of re¬ 
solved exoplanet being transited by an exomoon; initial 
spectroastrometric detection could be used to forecast 
the timing of these events. 

Spectroastrometry may also enable a measurement of 
the spectrum of the moon. In the case of an Earth-like 
planet orbiting a giant planet in the habitable-zone, this 
may be one of the best ways to characterize the atmo¬ 
sphere of the planetary companion. Once the solution 
for the moon-planet orbit is found, the observations can 
be integrated for a longer time, pos sibly allowing a spec¬ 
trum to be built up of both bodies ( p.2| ). This technique 
works best if two wave bands can be pinpointed within 
which the planet dominates in one band and the moon 
in the other; then the spectroastrometric signal between 
these two bands gives the sky separation of the bodies 
at each time. In t he e xample we considered of a face- 
on circular orbit (^4.2), the planet-moon separation is 
unchanged with time. Eor a different orbital configura¬ 
tion, the planet-moon separation needs to be modeled 
as a function of time and wavelength with a Keplerian 
orbit. Then the size of the orbit as a function of wave¬ 
length (with respect to a reference wavelength) can be 
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Figure 12. Minimum telescope diameters needed to obtain a 5—cr spectroastrometric detection between 0.45 fim and 0.89 fim of an 
Earth-like exomoon orbiting a warm Jupiter in the approximate habitable zones of the systems labeled above, displayed at their respective 
distances. The different colors represent the spectral types F (pale yellow), G (yellow), and K (orange). The circles represent the systems 
with confirmed exoplanets, the squares for systems with unconfirmed exoplanets, and the stars for systems in which exoplanets have yet 
to be detected. 


used for the spectral decomposition. The edge-on case 
will have a signal that is reduced by ^ \/2. This mea¬ 
surement can possibly allow the correction of a planet’s 
spectrum for contamination by a moon, avoiding a po¬ 


tential fals e-positive for disequilibrium chemistry (Rein 
et al.|[2M4 ). The recovered spectra of both bodies could 


constrain their albedos and their thermal emission prop¬ 
erties, further constraining the atmospheric properties. 

With measurements of the orbital and bulk properties 
of the moon-planet system, one can potentially place con¬ 
straints upon the formation of the system, as well as its 
tidal evolution. In particular, the location of the orbit 
within its Hill sphere, and the eccentricity of the orbit, 
will constrain tidal migration and circularization of the 
system, which can po tentially yield a co nstraint upon 
tidal evolution theory ( Heller et al. ||2Q14 ). 

We end by mentioning some other issues to consider 
in future work. We have only considered two exam¬ 
ple planet-moon compositions, and a limited range of 
an (admittedly) large parameter space. The sizes of 
the planet and moon, spectral type of the host star, 
planet-moon separation, spectral variability of the planet 
and/or moon, semi-major axis of the stellar orbit, exo- 
zodiacal light, presence of other planets, presence of rings 
(about the planet or star), inclination of the orbital axes, 
eccentricities, and more can affect the detection. Most 
importantly, the composition and presence (or absence) 
of an atmosphere, albedo, presence of clouds or hazes, 
and other properties of the atmospheres or reflective sur¬ 
faces will affect which bands should be observed, and 


at what wavelengths each body will dominate. Planet 
and/or moon variability should be considered as well. In 
the case of the Earth-Moon system, the quadrature Earth 
varies at the level of roughly ten percent over its rotation 
period (one day). This would be detectable in total flux, 
and so the spectrastrometric signal could be averaged 
over the rotation timescale, assuming that the two bands 
could be measured at the same time. Longer timescale 
variability could affect the orbital characterization and 
spectral disentangling; we leave the issue of variability 
for future work. Tidal heating of an exomoon could en¬ 
hance the spectroastrometric signal if the temperature 
were increased to the point of contributing strongly in 
the absorbing bands of the planet. Multiple moons (or 
rings) could affect the interpretation of the spectroastro¬ 
metric signal; we have assumed a single moon based on 
the logic that the first detections will be of large moons, 
and large moons might not exist stably with multiple 
moon companions. The inclination of the moon with re¬ 
spect to the observer will affect the spectroastrometric 
signal: an edge-on moon will spend some time at a small 
projected separation from the planet, and so an observa¬ 
tion at the wrong time might miss the spectroastromet¬ 
ric signal. This will require additional observing time 
to sample the (unknown) phase of the moon sufficiently 
to be able to capture it at maximum elongation. We 
have only considered observations at quadrature illumi¬ 
nation, while the orbital phase dependence of the spec¬ 
troastrometric signal may yield additional information 
about the bodies. Also, in computing the spectroastro- 
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metric signal, we have ignored the motion of the moon. 
This should be adequate if the duration of the observa¬ 
tion is much shorter than the orbital period of the moon, 
but for long exposure times or short moon periods, this 
could affect our signal-to-noise computations. Once a 
moon is detected, and its orbital period measured, the 
orbital motion can be accounted for in modeling the spec¬ 
troastrometric signal as a function of time. 

The design of the detector which would allow spec¬ 
troastrometric measurements needs to be explored. New 
detector technologies that allow for the measurement of 
the position and energy of a photon could potentially 
yield low-resolution spectroscopic images, and the result¬ 
ing data cube could be binned in wavelength channels 
to optimize the detection of a moon, and later to opti¬ 
mize the characterization of the spectra as a function of 
wavelength. Specifically, Microwave Kinetic Inductance 
Detectors, or MKIDs, and Superconducting Tunnelling 
Junctions, STJs, have high quantu m efficiency and are 
being multiplexed into larger arrays (Peacock et al.|19^ 


Day et al. 2003; Mazin et al. 2012). Another possible 
approach would be to use an integral field unit spec- 
trogr aph, such as that used on the Gemini Planet Im¬ 
ager (Macintosh et al. 2006). The wave bands that are 
used could in principle be selected after the observations 
are made: the smaller body tends to dominate at wave¬ 
lengths where the larger body is dimmest, so the wave¬ 
bands can be selected based on the observed molecular 
band features that are found after the initial spectrum is 
measured. Then the centroid offset can be searched for 
by measuring on-band and off-band. The spectroastro¬ 
metric precision may be affected by the sampling of the 
point-spread function as a function of wavelength; this 
should be studied further to determine what pixel/fiber 
size is required to obtain astrometric precision at each 
wavelength that is limited by the width of the PSF and 
number of photons detected, and not its spatial sampling. 

We have limited our exploration to a coronagraphic 
telescope operating out to 3 /im. Interferometric 
approaches might allow spectroastrometry to be ex¬ 
tended to longer wavelengths, where additional absorp¬ 
tion bands could be used and where the thermal emission 
from the moon can be significant, such as in the TPF 


I/Darwin concept operating from 6-20 jam (Lawson et 
^|2QQ8 Cockell et al.J|2QQ9| . This would have great ad¬ 


vantages for characterizing the planet-moon system, but 
involves the significant technical challenges of formation 
flying, baseline control, and thermal backgrounds. In 
principle this detection technique could be applied to 
ground-based direct imag ing as well with future large 
coronagraphic telescopes; Schneider et al. (2015) have 
mentioned the possibility of detecting the astrometric 
wobble of a directly imaged giant planet due to an ex¬ 
omoon companion. The precision may be limited more 
strongly by control of the optics and atmosphere, but 
nonetheless could possibly detect large moons in the in¬ 
frared, and spectroastrometry could potentially be more 
sensitive than absolute astrometry depending on the in¬ 
strumental and observational design. 

One challenge of coronagraphic spectroastrometry is 
that precise centroiding requires controlling sources of 
noise, as well as a stable PSF that can be measured astro- 
metrically over a broad range of wavelengths, requiring 
stable pointing and control of the PSF shape. These re¬ 


quirements may drive stricter tolerances for the design of 
a high-contrast imaging telescope than the requirements 
for spectroscopy. For example, the statistical uncertainty 
for the 28-day observation of a 1.34-pc Earth-Moon sys¬ 
tem with a 12-meter telescope is ~ 3 /iarcsec, and the 
typical ratio of the astrometric precision can be as small 
as 0.01-0.1% of the width of the PSF. Thus, the point¬ 
ing, the variation of the PSF, the characterization of the 
PSF, and the effects of other sources of noise need to be 
limited/controlled/calibrated to a very precise level. The 
telescope design and detector calibration would need to 
be capable of reaching this sub-pixel precision in measur¬ 
ing the location of the centroid of the planet-moon PSF 
as it passes across the detector over time; sub-pixel sen¬ 
sitivity variations, detector latency, charge bleeding, and 
light loss would all have to be controlled or calibrated 
to high precision to reach the photon-noise limit that we 
have assumed for this paper. One mitigating factor is 
that the spectroastrometric signal changes sign and di¬ 
rection as the moon orbits the planet, and should repeat 
with time, while any systematic errors will likely have a 
different time dependence. 

6. CONCLUSIONS 

Spectroastrometric detection presents a scientifically 
promising, but technically challenging, means of detect¬ 
ing and studying satellites of exoplanets, ‘exomoons.’ We 
have presented two case studies to illustrate the poten¬ 
tial for detection of this effect, and urge that studies of 
future coronagraphic telescope designs take into account 
the spectroastrometric signal when designing the techni¬ 
cal requirements of telescopes and detectors. Subsequent 
to detection, the characterization of exomoons with spec¬ 
troastrometry could lead to measurement of the funda¬ 
mental properties of exoplanet-moon systems, including 
the mass of the host planet, and could help in pinpointing 
targets that are valuable for studies of habitability. 

We end by making some general recommendations for 
guiding design of future coronagraphic space telescopes 
with spectroastrometry in mind: 1) the instrument suite 
should include the capability of making astrometric mea¬ 
surements as a function of wavelength for a range of 
spectral resolutions; 2) the pointing control, PSF wing 
suppression, and PSF stability and calibration should be 
designed with the capability of making photon-limited 
astrometric measurements over a broad range of wave¬ 
lengths; and 3) the instrument wavelength should be ex¬ 
tended out to ~ 3 micron to cover the water absorption 
band at which the Moon dominates over the Earth (« 2.7 
micron) and methane band at which the Earth dominates 
over Jupiter (^2.3 micron). 
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Table 2 

Highest SNR combination of wavelengths for a representative selection of spectral resolutions for the Earth-Moon system set at 1.34 
parsecs from the observer using a 12-m coronagraphic space telescope exposed for 24 hours. 


R 

(A/AA) 

SNR 

Signal 

(mas) 

Ai 

(fim) 

AAi 

{fim) 

Moon Flux 
% of Total 

Planet Flux 
% of Total 

A 2 

(fim) 

AA 2 

(/rm) 

Moon Flux 
% of Total 

Planet Flux 
% of Total 

1.0 

17.86 

0.1264 

0.295 

0.147-0.442 

0.75 

99.25 

1.998 

0.999-2.997 

7.34 

92.66 

1.5 

19.30 

0.2138 

0.328 

0.219-0.438 

0.74 

99.26 

2.010 

1.340-2.681 

11.89 

88.11 

2.0 

18.56 

0.4082 

0.344 
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99.28 

2.390 

1.792-2.987 
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9.43 
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4.0 
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0.60 
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5.0 
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0.348 
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1.939 

1.920-1.958 

58.29 

41.71 

60.0 

7.86 

1.1152 

0.349 

0.346-0.352 

0.49 

99.51 

1.947 

1.931-1.963 

58.64 

41.36 

100.0 

6.35 

1.1622 

0.347 

0.346-0.349 

0.46 

99.54 

1.949 

1.940-1.959 

61.07 

38.93 

150.0 

5.28 

1.1572 

0.349 

0.348-0.350 

0.45 

99.55 

1.955 

1.949-1.962 

60.80 

39.20 

200.0 

4.73 

1.1759 

0.348 

0.347-0.349 

0.57 

99.43 

1.931 

1.926-1.936 

61.89 

38.11 


Table 3 

Highest SNR combination of wavelengths for a representative selection of spectral resolutions for the Jovian-Earth system set at 10 
parsecs from the observer using a 12-m coronagraphic space telescope exposed for 24 hours. 


R 

(A/AA) 

SNR 

{^) 

Signal 

(mas) 

Ai 

(fim) 

AAi 

(fim) 

Moon Flux 
% of Total 

Planet Flux 
% of Total 

A 2 

{fim) 

AA 2 

{fim) 

Moon Flux 
% of Total 

Planet Flux 
% of Total 

1.0 

3.33 

0.0088 

0.858 

0.429-1.286 

0.54 

99.46 

1.288 

0.644-1.931 

0.97 

99.03 

2.0 

11.99 

0.0450 

0.571 

0.428-0.714 

0.33 

99.67 

0.947 

0.710-1.184 

2.53 

97.47 

3.0 

15.07 

0.0775 

0.514 

0.428-0.600 

0.25 

99.75 

0.858 

0.715-1.001 

4.03 

95.97 

4.0 

16.59 

0.1282 

0.490 

0.429-0.552 

0.21 

99.79 

0.887 

0.777-0.998 

6.47 

93.53 

5.0 

16.53 

0.1680 

0.476 

0.428-0.523 

0.20 

99.80 

0.899 

0.809-0.989 

8.40 

91.60 

10.0 

16.01 

0.3812 

0.452 

0.429-0.474 

0.18 

99.82 

0.925 

0.879-0.971 

18.79 

81.21 

20.0 

12.64 

0.3119 

0.456 

0.445-0.468 

0.17 

99.83 

0.882 

0.860-0.904 

15.40 

84.60 

30.0 

13.49 

0.5779 

0.451 

0.443-0.458 

0.17 

99.83 

0.894 

0.879-0.908 

28.38 

71.62 

40.0 

15.55 

1.0089 

0.449 

0.443-0.454 

0.16 

99.84 

0.892 

0.881-0.903 

49.42 

50.58 

50.0 

15.64 

1.2301 

0.452 

0.447-0.456 

0.16 

99.84 

0.891 

0.882-0.900 

60.22 

39.78 

60.0 

14.73 

1.2646 

0.455 

0.451-0.459 

0.16 

99.84 

0.890 

0.882-0.897 

61.91 

38.09 

80.0 

13.25 

1.3377 

0.448 

0.445-0.451 

0.15 

99.85 

0.890 

0.884-0.896 

65.47 

34.53 

100.0 

12.20 

1.3499 

0.447 

0.445-0.450 

0.15 

99.85 

0.887 

0.883-0.892 

66.06 

33.94 

120.0 

11.04 

1.3296 

0.453 

0.451-0.455 

0.14 

99.86 

0.887 

0.883-0.891 

65.06 

34.94 

150.0 

9.86 

1.3327 

0.470 

0.468-0.471 

0.13 

99.87 

0.890 

0.887-0.893 

65.21 

34.79 

200.0 

8.53 

1.3319 

0.430 

0.429-0.431 

0.12 

99.88 

0.886 

0.884-0.889 

65.15 

34.85 

250.0 

7.68 

1.3217 

0.491 

0.490-0.492 

0.12 

99.88 

0.889 

0.887-0.891 

64.66 

35.34 

300.0 

7.22 

1.3790 

0.431 

0.430-0.432 

0.10 

99.90 

0.889 

0.888-0.890 

67.43 

32.57 

400.0 

6.32 

1.4111 

0.431 

0.430-0.431 

0.09 

99.91 

0.890 

0.889-0.891 

68.99 

31.01 



